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Abstract 

We calculate the relation between the bare coupling constant go and the MS- 
renormalized coupling ffj^g, go = Z g (go, ap)g^, to 2 loops in perturbation 
theory. We employ the standard Wilson action for gluons and the overlap 
action for fermions. For convenience, we have worked with the background 
field technique, which only requires evaluation of 2-point Green's function 
for the problem at hand. Our results depend explicitly on the number of 
fermion flavors (Nf) and colors (N). Since the dependence of Z g on the 
overlap parameter p cannot be extracted analytically, we tabulate our results 
for different values in the allowed range of p (0 < p < 2), focusing on values 
which are being used most frequently in simulations. Knowledge of Z g allows 
us to derive the 3-loop coefficient of the bare /3-function (/3l(<?o)) which, unlike 
the 1- and 2-loop coefficients, is regularization-dependent. The nontrivial 
dependence of Z g and of /3l(<?o) on p is plotted for various choices of N, Nf. 

Keywords: Beta function, Overlap fermions, Running coupling constant, 
Lattice perturbation theory. 
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I. INTRODUCTION 



In later years, use of non-ultralocal actions which preserve chiral symmetry on the lattice 
has become more viable. The two actions which are being used most frequently are overlap 
fermions [1-3] based on the Wilson fermion action and domain-wall fermions [4,5]. 

Overlap fermions are notoriously difficult to study, both numerically and analytically. 
Many recent promising investigations involving simulations with overlap fermions have ap- 
peared; see, e.g., Refs. [6-12]. Regarding analytical computations, the only ones performed 
thus far have been either up to 1 loop, such as Refs. [13-19], or vacuum diagrams at higher 
loops [20,21]. The present work is the first one involving non-vacuum diagrams beyond the 
1-loop level. 

We compute the 2-loop renormalization Z g of the bare lattice coupling constant go in 
the presence of overlap fermions. We relate go to the renormalized coupling constant gy[§ 
as defined in the MS scheme at a scale £l; at large momenta, these quantities are related as 
follows 

a Ms(f L ) = + di(p,a)al + d 2 (fia)a i + ... , (1) 

(a = gll^K-, «ms = fi'Ms/^ 71 "' a '■ lattice spacing). The 1-loop coefficient di(fla) has been 
known for a long time; several evaluations of <i 2 (/ia) have also appeared in the past ~10 
years, either in the absence of fermions [22,23], or using the Wilson [24] or clover [25,26] 
fermionic actions. Knowledge of c^Om), along with the 3-loop MS-renormalized /3-function 
[27] allows us to derive the 3-loop bare lattice /3-function, which dictates the dependence of 
lattice spacing on go. In particular, it provides a correction to the standard 2-loop asymptotic 
scaling formula defining A^. Ongoing efforts to estimate the running coupling from the lattice 
[28-31] have relied on a mixture of perturbative and non-perturbative investigations. As a 
particular example, relating to a$F (SF: Schrodinger Functional scheme, as advocated 
by the ALPHA Collaboration), entails an intermediate passage through the bare coupling 
and the conversion from a^g to «o is carried out perturbatively. 

The paper is organized as follows: Some theoretical background and the methodology 
of the necessary perturbative calculations are given in Section II. Section III regards the 
overlap action where the derivation of the vertices with up to 4 gluons is explained. We 
also provide the expressions for these vertices (the 4-gluon vertex is given in Appendix A). 
Details on our computation, numerical results and plots of Z g coefficients and the /9-function 
can be found in Section IV. Finally, in Section V we give the summary and our conclusions. 

The present study, being the first of its kind in calculating 2-loop diagrams with overlap 
vertices and external momentum dependence, had a number of obstacles to overcome. One 
first complication is the size of the algebraic form of the Feynman vertices; as an example, 
the vertex with 4 gluons and a fermion- ant if ermion pair contains ~724,000 terms when ex- 
panded. Upon contraction these vertices lead to huge expressions (many millions of terms); 
this places severe requirements both on the efficiency of the computer algorithms which we 
must design to manipulate such expressions automatically and on the necessary computer 
RAM. Numerical integration of Feynman diagrams over loop momentum variables is per- 
formed on a range of lattices, with finite size L, and subsequent extrapolation to L — > oo. 
As it turned out, larger L were required for an accurate extrapolation in the present case, 
compared to ultra-local actions. In addition, since the results depend nontrivially on the 
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parameter p of the overlap action, numerical evaluation must be performed for a sufficiently 
wide set of values of p, with an almost proportionate increase in CPU time. Extreme values 
°f P (P ^ 0) P ^ 2) show unstable numerical behaviour, which is attributable to the spu- 
rious poles of the fermion propagator at these choices; this forces us to even larger L. A 
consequence of all complications noted above is an extended use of CPU time: Our numer- 
ical integration codes, which ran on a 32-node cluster of dual CPU Pentium IV processors, 
required a total of ~50 years of CPU time. 



II. THEORETICAL BACKGROUND 



The definition and value of the renormalized coupling constant g depends on the renor- 
malization scheme (parameterized by a scale p) , and this dependence is given by the renor- 
malized /5-function 



P(g) = p 



dg_ 

dp 



(2) 



We will adopt the MS renormalization scheme throughout this work, and we will denote 
the renormalized coupling constant simply by g. For the lattice regularization a bare 
/3-function is defined as 



(3) 



where p the renormalization scale and g (go) the renormalized (bare) coupling constant. It 
is well known that in the asymptotic limit for QCD (g — > 0), one can write the expansion 
of the /5-function in powers of g , that is 

Pl (9o) = ~bo 9o ~ h g 5 -b^g 7 -... (4) 
The coefficients b , b\ are universal constants (regularization independent) given by 

1 
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N — N 



f 



3 N 



(5) 



where N(Nf) is the number of colors (flavors). 

On the contrary, bf (i > 2) depends on the regulator; it must be determined perturbatively. 
In the present work we calculate the coefficient b\ using the overlap fermionic action and 
Wilson gluons. 

The perturbative expansion of the renormalized /9-function is 



«»>- 4 



-bo g 3 -h g 5 -b 2 g 7 + .. 



(6) 



a,9o 
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Al(<?o) and (3(g) can be related using the renormalization function Z g , defined through 
9o = Zg(g ,ajj,)g, that is 1 



1 " % —Q^r- J Z g P(9oZ g ) 



Computing Z 2 to 2 loops 



Z g 2 (^ , a/x) = f + So 2 (26 Haji) + Z ) + <7o (2&i ln(a/i) + Zi) + O(^ 6 ) 



(7) 



(8) 



and inserting it in Eq. (7), allows us to extract the 3-loop coefficient b 2 . The quantities 
bo,bi,b 2 and Iq have been known in the literature for quite some time [27,19]; &o and b\ are 
the same as those of the bare /3-function, Eq. (5), and b 2 in the MS scheme is 
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The constant Iq is related to the ratio of the A parameters associated with the particular 
lattice regularization and the MS renormalization scheme 



l = 2b \n (Al/Ams) 
For overlap fermions the exact form of Iq appears in Ref. [19] 



l = — - 0.16995599iV + N f 
8N 1 



72tt 2 



-k(p) 



(10) 



(11) 



where k(p) is the convergent part of the 1-loop fermionic contribution (denoted by in 
Ref. [19]), presented in Table I. 

Eq. (7) is valid order by order in perturbation theory and expanding it in powers of g\ 
the first nontrivial relation is 



b 2 - Mo + hh 



(12) 



Thus, the evaluation of b^ requires only the determination of the 2- loop quantity l ± . 

A direct outcome of our calculation is the 2-loop corrected asymptotic scaling relation 
between a and g 



A, 



exp 



2&o<?6\ 



{bQgl)- h ^\l+qgl + 0{gt) 



Q 



b\ - bob 



2bl 



(13) 



where all quantities in the correction term q, except b 2 , are known. 

The most convenient and economical way to proceed with calculating Z g (g , ajj) is to 
use the background field technique [32,33,23], in which the following relation is valid 



1 Z g could be denoted: Z g ' MS to indicate its dependence on the regulator (L: lattice) and on the 
renormalization scheme (MS). 
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Za{9o, a/i)Zg(go, a/i) = 1 (14) 

where Za, defined as: A^(x) = ZA(go, ap') 1 ^ 2 A R (x) is the background field renormalization 
function (A^ (A R ): bare (renormalized) background field). In the lattice version of the 
background field technique, the link variable takes the form 

U^x) = e ia9oQ ^ {x) ■ e iaA » {x) (15) 

(Q^x): quantum field, A^x): background field). In this framework, instead of Z g (gQ,ajl), 
it suffices to compute ZA(go,ap,), with no need to evaluate any 3-point functions. For this 
purpose, we consider the background field one-particle irreducible (1PI) 2-point function, 
both in the continuum (dimensional regularization, MS subtraction): r^ A (p)"^ and on the 
lattice: T AA (p)^ u . In the notation of Ref. [23], these 2-point functions can be expressed in 
terms of scalar functions vr(p), v{p) 

r* 1 ^ = (W 2 - Mu) (i - Mp)) M Mp) = 9 2 4\p) + 9 4 4\p) + - (16) 

Y.Tt A {p)i = S ab 3p 2 [1 - u(p)] /gl u(p) = gl^ip) + <? V 2 )(p) + ... (17) 

(ftp = (2/a) sin(ap M /2)). There follows 

1 - unfa ft, g) 

^ A = 1 7 7 ^ 

1 - v(p, a, g ) 

The gauge parameter A must also be renormalized (up to 1 loop), in order to compare lattice 
and continuum results 

X = Z Q X , Z Q = l+g 2 z$ ) + ... (19) 

{Zq\ renormalization function of the quantum field). The coefficient Zq 1 is obtained from 
the quantum field 1PI 2-point function in the continuum (r^ Q (p)^) and on the lattice 
^ Q L Q {p)t) through 



r? Q (p)J = S ab [(<W 2 - P&) (1 " MP)) + ^P,Pu] , Mp) = 9 2 " { r\p) + 0(g 4 ) (20) 
E r L Q G< = -* ab P 2 [3 (1 - v(p)) + Ao] , u(p) = 9l^\p) + 0{gt) (21) 



4 1} =uW(p,a,g )-u ( £>( P ,il,g) (22) 
In terms of the perturbative expansions Eqs. (16), (17), (20), (21), Z 2 g takes the form 

A=A 



a (1) 

1 + 9l (4 ] - + 9t (4 2) - - (2) ) + W " 4 ] )-^} (23) 



(i) 



The fermion part of coincides with that of Similarly for the fermion part of uj r 
and Vr 1 . Consequently, one may write 
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(24) 



Since the quantities of interest are gauge invariant, we choose to work in the bare Feyn- 
man gauge, A = 1, for convenience. In order to compute Za we need the expressions for 
v r\ u r\ z q-> u ^-> The MS renormalized functions necessary for this calculation to 2 
loops are [23,34] 
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i/£>(p,A=l) = 
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(27) 



For the lattice quantities, the gluonic contributions (Nf = 0) have been presented in 
previous works [23,34] (for the Wilson action) 



CO 



(1) 



V 



(1) 



(p,A„ = l) 



(P,A = 1) = 



In (aV) - + 0.137286278291N 



^ (2) (P ; A = 1) 



lliV 
'48^ 

N 
"32tt 4 



In (a V) 



ln(aV) + 



8N 
3 



128N 2 



+ 0.217098494367N 



- 0.01654461954 + 0.0074438722^ 



(28) 



(29) 



(30) 



The fermionic contributions are associated with the diagrams of Fig. 1 and Fig. 2. In the 
present work, is perturbatively calculated for the first time using overlap fermions and 
Wilson gluons. For completeness, we also compute the coefficient and compare it with 
previous results. The 1-loop diagrams (Fig. 1) correspond to u^, while the 2-loop diagrams 
(Fig. 2) lead to u^ 2 \ Dashed lines ending in a cross represent the background field, while 
those inside loops denote the quantum field. Solid lines correspond to fermion fields and a 
dot stands for the mass counterterm. Note that, for overlap fermions, the mass counterterm 
equals zero, by virtue of the exact chiral symmetry of the overlap action; consequently, 
diagrams 19 and 20 both vanish. Certain 2-loop diagrams have infrared divergences and 
become convergent only when grouped together (6+12, 7+11, 8+18, 9+17). 
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Fig. 1: Fcrmion contributions to the 1-loop function i/ 1 ). 
ground gluons. Solid lines represent fermions. 



Dashed lines ending on a cross represent back- 




5 6 7 8 




9 10 11 12 




13 14 15 16 




17 18 19 20 

Fig. 2: Fermion contributions to the 2-loop function v^ 2 \ Dashed lines represent gluonic fields; those ending 
on a cross stand for background gluons. Solid lines represent fermions. The filled circle is a 1-loop fermion 
mass counterterm. 
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III. OVERLAP ACTION 



In recent years, overlap fermions are being used ever more extensively in numerical simu- 
lations, both in the quenched approximation and beyond. This fact, along with the desirable 
properties of the overlap action, was our motivation to calculate the /9-function with this 
type of fermions. The important advantage of the overlap action is that it preserves chiral 
symmetry while avoiding fermion doubling. It is also 0(a) improved. The main drawback 
of this action is that it is necessarily non-ultralocal; as a consequence, both numerical simu- 
lations and perturbative studies are extremely difficult and demanding (in terms of human, 
as well as computer time). 

The overlap action is given by [2] 



^overlap = CL % ^ D N (n, TJl) *(m) 



(31) 



where D N (n,m) is the overlap-Dirac operator 



D N (n,m) = p 



- X 



Vxvc „ 



X 



[Dw - p) 



(32) 



and Dw is the Wilson-Dirac operator 
1 



7„ (V; + V„) - aV^Vj , V^{x) = - [U(x, p)^(x + ap) - ^{x)\ (33) 



The overlap parameter p is restricted by the condition < p < 2 to guarantee the correct 
pole structure of D N . The coupling constant is included in the link variables, present in the 
definition of X, and one must take the perturbative expansion of X in powers of g . This 
expansion in momentum space takes the form 

X(p',p) = X o(p)(2n) 4 5 P (p' - p)+X 1 (p',p) + X 2 (p',p) +X 3 (p',p) + X 4 (p',p) +O(g 5 ) (34) 



tree— level 



1—loop 



2— loop 



where Xo is the inverse fermion propagator and Xi are the vertices of the Wilson fermion 
action (p (p'): fermion (antifermion) momentum). The construction of all overlap vertices 
relevant to the present computation (see Eqs. (39), (40), (41), (42), (A7) below) make use 
of Xo and X\ — X 4 ; these quantities can be written in the compact form 



Xoip) = - H 7 M sin(ap M ) + - ^ (l - cos(ap M )) - - 
X 1 (p',p)=g J d'kSip'-p-k^A^V^^ 



d%d4h 5(p'-p- ki - k 2 )J2AMAMV 2 „(^) 



Xs(p',p) 



9l 



9j 
4! 



d 4 kid 4 k 2 d 4 k 3 



ti i=i z 



i=i 



d 4 kxd 4 k 2 d 4 k 3 d 4 k 4 



5( P > - P - e e n Mb) [ - « 2 M^)1 ( 35 ) 



(2*) 
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where 

Vi,fj,(p) = z7 M cos(ap M ) +r sir^ap^), V2 )At (p) = -27^0 sh^ap,,) + ar cos(ap M ) (36) 

A M represents a gluon field; later on we will have to generalize Eqs. (35) to the case where 
both background and quantum gluon field are present, see Eq. (43). 

At this point we can proceed with the perturbative expansion of D N in powers of g . 
This leads to the propagator of zero mass fermions and to gluon-fermion-antifermion vertices 
(with up to 4 gluons for the needs of the present work). The much simpler case of vertices 
with up to 2 gluons (and no background) can be found in Ref . [35] . After laborious analytical 
manipulations (an essential step is the expansion of i/Vxvc using complex analysis, which 
is presented in Appendix A), the overlap-Dirac operator is expanded into terms with up to 
4 gluons as 

D N (k x , k 2 ) = Dofa) (2vr) 4 5\k x - k 2 ) + E(/ci, k 2 ) (37) 
D (ki) is the inverse propagator, 



Do(h) = 1 + w(p) = ./(£sin 2 (p M )) + (p ~ 2r £sin 2 (^/2)) Z (38) 



and 



^hlM= yifa fc 2 ) + y2 {ku h) + V 2 {ku fca ) 

1— gluon vertex 2— gluon vertex 

+ v?(k h k 2 ) + v£(h, k 2 ) + v£(h, k 2 ) 

v ' 

3 — gluon vertex 

+ Vi(ki, k 2 ) + V 2 \h, k 2 ) + V£(h, k 2 ) + Vi(k h k 2 ) +0{gl) (39) 

4— gluon vertex 

where we have set a = 1. Vl,V 2 % ^ ^3 are given below and the reader can find the expression 
for V£ in Appendix A. 



vi{h, k 2 ) 



u(ki) + u;(k 2 ) 



Xi(h, k 2 ) 1 r Xo(fci) X}(k 1} k 2 ) xo(h 

oj(ki)u)(k 2 ) 



(40) 



- r d 4 k 3 1 1 1 

2( u 2) ~ J (2^p ufa) + u(k 3 ) u)(ki) + u(k 2 ) uj{k 2 ) + 



u;(/c 3 



x 



J2 - Xj(ki, k 3 ) xl(h) X k (k 3 , k 2 ) 

{j>0,fe>0} L 

{j+k=i} -Xjih, k 3 ) Xl(k 3 , k 2 ) X o{k 2 ) 
— Xo(fci) Xj(kt, k 3 ) X k (k 3 , k 2 ) 

+ ^(fe 1 ) + C,(fc 2 ) + ^(fc3) x t (fci; h)Uh)X l {k ^ k2)Mk2 

uj{ki)uj{k 2 )uj{k 3 ) 



(41) 
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d 4 k 3 d 4 k^ 1 



n 



{j>o,fc>o,i>o} L 

{j+k+l=i} 



(27r) 4 (2n) 4 A\ p ^ 4 uj(k Pl )+uj(k P2 )J 

u ( k Pi) u ( k P2) UJ (k p . i UX j (k 1 , h)Xl(k 3 ,k 4 )Xi(k 4 ,k 2 ) 

" \peS4 / 



+ (co(k 1 ) +u(k 3 ) +u(k 4 ) +co(k 2 )) x 

Xo(*i)-Xj(*i, h)X k (k 3 , fc 4 )X/(fc 4 , k 2 )xo(k 2 ) + Xo{h)X](ki, k 3 )X k (k 3 , k 4 )xl (tyX^k*, h) 
+Xo(h)Xj(k 1: k 3 ) Xo (k 3 )Xl(k 3 , k^X^h, k 2 ) + X^/cx, k 3 )xl(k 3 )X k (k 3: h)Xj(k 4} k 2 ) Xo (k 2 ) 
+X j (k 1 , k 3 )xl(k 3 )X k (k 3 , k A )xl{ki)Xi{ki, k 2 ) + X^, k 3 )Xl(k 3 , h) X o(h)Xj(h, k 2 ) X o(k 2 ) 
u(k Pl )uj(k P2 )(uj(k Pl )/2 + u;(A;p 3 )/3 



(42) 



V P 5 4 uj(h)uj(k 3 )uj(h)uj(k 2 ) J " 

Xo(A;i)Xj(/c 1 , A; 3 )xo(A;3)^i(^3, fc4)Xo(k4)^/(&4, k 2 )xo(k 2 ) 



(£4: permutation group of the numbers {1,2,3,4}) 

The use of the background field technique implies that instead of the generic gluonic fields 
(appearing in X{S, Eq. (35)), one must consider all possible combinations of background 
(A) and quantum (Q) fields which originate in the links of Eqs. (33), (15). Hence, 

X 1 {p',p) = X?{p',p)+X A {p',p) 

X 2 (p',p) = X? Q (p',p) + Xf \p',p) + Xi A (p',p) 

Xs(p',p) = X? QQ (p>,p) + X$ QA {p',p) + X^ AA (p',p) + X AAA (p',p) 

X,(p',p) = X? QQQ (p>,p) + X? QQA (p>,p) + X? QAA (p',p) + Xf AAA (p',p) + X AAAA (p',p) (43) 
As an example, let us write the expression for X® QA 



X^ QA (^' 



(p',p) '- 

E 



gl f d 4 k\d 4 k 2 d 4 ki 

'p' + p 
~2 



5{p' -p-ki-k 2 -k 3 ) x 



a 2 V 1 ^(^^)(Q, l (k 1 )Q^(k 2 )A^k 3 ) + A^k 3 )Q, L (k 2 )Q p (k, 



(44) 



Upon substituting the expression for Xj's in the overlap vertices, the latter become extremely 
lengthy and complicated. For instance, the vertex with Q-Q-A-\l/-\l/ consists of 9,784 terms, 
while the vertex with Q-Q-A-A-\I>-* has 724,120 terms. 



IV. RESULTS 

For the algebra involving lattice quantities, we make use of our symbolic manipulation 
package in Mathematica, with the inclusion of the additional overlap vertices. The first step 
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to evaluate the diagrams is the contraction among vertices, a step performed automatically 
once the vertices and the 'incidence matrix' of the diagram are specified. The outcome of 
contraction is a preliminary expression for the diagram under study; there follow simplifica- 
tions of the color dependence, Dirac matrices and tensor structures. We use symmetries of 
the theory (permutation symmetry and lattice rotational invariance), or any other additional 
symmetry that may appear in particular diagrams, to keep the size of the expression down 
to a minimum. The external momentum p appears in arguments of trigonometric functions 
and the extraction of p dependence is divided into two parts: first we isolate terms that 
give single and double logarithms (a few thousand terms, expressible in terms of known, 
tabulated integrals), and then for the convergent terms we employ naive Taylor expansion 
with respect to p up to 0(p 2 ). This extraction makes explicit the functional dependence of 
each diagram on p; the coefficients of terms proportional to p 2 are integrals over the two 
internal momentum 4-vectors. The required numerical integrations are performed by opti- 
mized Fortran programs which are generated by our Mathematica 'integrator' routine. Each 
integral is expressed as a sum over the discrete Brillouin zone of finite lattices, with varying 
size L, and evaluated for different values of the overlap parameter p. The average length of 
the expression for each diagram, after simplifications, is about 2-3 hundred thousand terms, 
so that diagrams must be split into parts (usually of 2000 terms) to be integrated. The nu- 
merical values of these parts must then be added together to avoid running into systematic 
errors or spurious divergences. Finally, we extrapolate the results to L — > oo; this procedure 
introduces an inherent systematic error, which we can estimate quite accurately. Infrared 
divergent diagrams must be summed up before performing the extrapolation. 

We denote the contribution of the i th 1-loop Feynman diagram to v^(jp) as v^\p)'i 
similarly, contributions of 2- loop diagrams to v^ 2 \p) are indicated by v^\p). The quantities 
vpijp), vf\p) depend on N, Nf, p and ap according to the following formulae (A = 1) 



ap 2 vV(p) = Nf 



ap 2 vf\p) = N } 



kf+a 2 p 2 {k^ + kP^} + 0((apy) 



(4tt) 5 



cf ] + a 2 p 2 t 



(1U (2) lnaV (3) (\na 2 p 

c + cy - Nn + c) — — - + c. 



(Air) 



(4tt)< 



(4) H 



(p 



2\2 



+ 0{(apf) 



(45) 



(46) 



where p 2 = 4 sin 2 (p M /2). The index % runs over diagrams, and the coefficients k\ , c 



depend on the overlap parameter p. Moreover, c~P = \cf' L> /N + cf' 1 ' N . Comparison with 
continuum results and usage of Ward Identities requires 

• ^ = ^5 X] c «'°' ) = (gauge invariance) 



c- 4 '* = (Lorentz invariance) 



(3) 1 (3) _ 4 AT (3) 

# C 15 _ gjy) C 16 _ g iV ) C 17 



(3) 



N 2 -l 
3N 



(47) 
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We have checked that all the above conditions are verified by our results. Inserting these 
conditions in Eqs. (45), (46), the expressions for the fermionic contribution to u^(p), u^(p) 
(after addition of all diagrams) take the form 

2 In a 2 p 2 



v 



(i) 



(p) 



(p) 



A) 



,( 2 ) 



(p) 



(p) 



7V /= 



3 (4tt) ; 



+ ^(M 2 ) 



N f =0 



.(i) 



(48) 



(49) 



In Table I we tabulate the total 1-loop contribution k^ = Y^ikY' for 21 values of the 
overlap parameter (0 < p < 2). Each diagram was integrated for lattice size L 4 , L < 128 
and then extrapolated to L — > oo. In all Tables and Figures, the errors accompanying 
our results are entirely due to this extrapolation. The coefficients k^\p) do not depend 
on the number of colors N, nor on the choice of regularization for the pure gluonic part 
of the action (Symanzik, Iwasaki, etc.). The numbers in Table I are in agreement with 
corresponding numbers from Ref. [19]. 

The 2-loop calculation of is^(p) was accomplished for the same 21 values of p and for 
L < 28. Table II presents the coefficients c^~^ = X)j cf 1 ' -1 ^, c^ 1 ' 1 ) = X)j c^ 1 ' 1 ^, versus p. Due 
to the extremely large size of the vertices involved, it is almost impossible to extend the 
results to larger L. Typically, the integration of 2000 terms is completed in ~ 7 days on 
1 CPU; the present calculation comprises approximately 3500x2000 terms. Thus, if only a 
single CPU were available, our work would have required ~50 years. In certain cases with 
large systematic errors we extended the integration up to L = 46. In general, the overlap 
action leads to coefficients which are very small for most values of p. As a consequence, 
systematic errors, which are by and large rather small, tend to be significant fractions of the 
signal for p > 1.4. 

From Eqs. (8), (23)-(28), (30) we find the following expression for l x in terms of N, Nf, 



k 



(l) (i,i) 



128 N 2 



+ 0.018127763034 - 0.007910118514 N 2 



+N f 



(167r 2 ) 2 iVV12 



(5|-4C(3)) 



N 481 



N 



(16tt 2 ) 2 36 8tt 2 



-(i-i) 



(50) 



We can write the final form of the 3-loop coefficient b\ for the /9-function (Eq. (12)), including 
gluonic as well as fermionic contributions, using Eqs. (5), (9), (11) and (50) 
11 



b L 2 = 



2048tt 2 N 



+ 0.000364106020 N - 0.000092990690 N a 



+N f 



{An 2 - l) 2 
4(167T 2 ) 3 A^ 2 

Nf ( 



0.000046883436 - 000013419574 N 2 



23 



(167T 2 ) 3 

(11 N - 



9N 3N 



8C(3) + 37 N" 



6 



2N f 



48tt 2 



.(i,-i) 



+ c 



(i,i) 



N) + 



(4N 3 + N f -3 N 2 N f 

(167T 2 ) 2 AT 



(51) 
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A large variety of possible numerical checks has been performed, as mentioned above: a. 
The total contribution to the gluon mass adds to zero, as expected, b. The coefficients of 
the non-Lorentz invariant terms cancel, c. The terms with double logarithms correspond 
to the continuum counterparts. This has been checked diagram by diagram, d. Terms with 
single logarithms add up to their expected value, which is independent of p (although the 
expressions per diagram are p-dependent). 

In Fig. 3 we plot the 1-loop coefficient k^ 1 ' with respect to p. Note that the errors are too 
small to be visible at this scale. The 2-loop coefficients c 1 - 1 ' -1 ^ and c^ 1 ' 1 ' are plotted in Figs. 
4-5, respectively, for different values of the overlap parameter. The extrapolation errors are 
visible for p < 0.4 and p > 1.7. Substituting fcW, c^ 1 ' -1 ) and c^- 1 ' 1 ' into Eq. (51), we find 
the numerical results for the 3-loop contribution, b\ , of the /3-function. These are plotted 
in Fig. 6, choosing N = 3 and Nf = 0, 2, 3. 




0.006 



Fig. 3: Plot of the total 1-loop coefficient = fcj 1 ' versus the overlap parameter p 
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-0.008 



-0.010 



-0.012 
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P 

Fig. 4: Plot of the total 2-loop coefficient c^ 1 ' -1 ^ = c ^~ 1 " 1 versus p. 



0.005 




Fig. 5: Plot of the total 2-loop coefficient c^ 1 ' 1 ^ = ^ c- 1 versus p. 
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0.002 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

P 

Fig. 6: The 3-loop coefficient b\ (Eq. (51)), plotted against p, for N = 3 and Nf — (horizontal red line), 
Nf = 2 (green line) and Nf = 3 (blue line). 

V. SUMMARY AND CONCLUSIONS 

We have calculated the 2- loop coefficient of the coupling renormalization function Z g , 
for the Yang-Mills theory with gauge group SU(N) and Nf species of overlap fermions. We 
used the background field method to simplify the computation; in this method there is no 
need of evaluating any 3-point functions. This is the first 2-loop calculation using overlap 
fermions with external momenta, and it proved to be extremely demanding in human and 
CPU time; this is due to the fact that we had to manipulate very large expressions (millions 
of terms) in intermediate stages. 

We used our numerical results of Z g to determine the 3-loop coefficient b\ of the bare 
lattice /3-function; the latter dictates the asymptotic dependence between the bare coupling 
constant g Q and the lattice spacing a, required to maintain the renormalized coupling at a 
given scale fixed. Knowledge of b% provides the correction term to the standard asymptotic 
scaling relation between a and g , via Eq. (13). 

The dependence of Z g and b\ on iV and Nf is shown explicitly in our expressions. On 
the other hand, dependence on the overlap parameter p cannot certainly be given in closed 
form; instead, we present our results for a large set of values of p in its allowed range. 

The 3-loop correction is seen to be rather small: This indicates that the perturbative 
series is very well behaved in this case, despite the fact that it is only asymptotic in nature. 
Furthermore, around the values of p which are most often used in simulations (1 < p < 1.6), 
fermions bring about only slight corrections to the 3-loop /3-function, even compared to pure 
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gluonic contributions, as can be seen from Fig. 6. 

The only source of numerical error in our results has its origin in an extrapolation to infi- 
nite lattice size. Compatibly with the severe CPU constraints, numerical 2-loop integration 
had to be performed on lattices typically as large as 28 4 , or even up to 46 4 in cases where an 
improved extrapolation was called for. An intermediate range for p (0.6 < p < 1.3) showed 
the most stable extrapolation error, and this may be a sign of their suitability for numerical 
simulations. 

As a by product of the present work, we have produced the lengthy expressions cor- 
responding to all overlap vertices which can arise in a 2-loop computation, and presented 
them in a rather compact form. Further computations of similar complexity, for example 
the 2-loop renormalization of operators in the overlap action (such as fermion currents), only 
require the vertices which we have presented here. 

Acknowledgements: This work is supported in part by the Research Promotion Foun- 
dation of Cyprus (Proposal Nr: ENIEX/0505/45). 
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TABLES 



p 


1.(1) 


0.1 


0.020377(7) 


0.2 


0.01581702(2) 


0.3 


0.0133504717(2) 


0.4 


0.0116910952(1) 


0.5 


0.0104621922(2) 


0.6 


0.0095058191(2) 


0.7 


0.00874441051(7) 


0.8 


0.00813753230(4) 


0.9 


0.00766516396(3) 


1.0 


0.00732057894(3) 


1.1 


0.00710750173(2) 


1.2 


0.00703970232(7) 


1.3 


0.0071425543(2) 


1.4 


0.0074569183(2) 


1.5 


0.0080467046(1) 


1.6 


0.0090134204(1) 


1.7 


0.010526080(2) 


1.8 


0.0128914(2) 


1.9 


0.01680(8) 



TABLE I. Coefficient k^ = ^ for different values of the overlap parameter 

i 



p 


c (i,-i) 


c (i,i) 


0.1 


-0.0096(6) 


0.124(3) 


0.2 


-0.0044(1) 


0.0118(5) 


0.3 


-0.00321(6) 


0.0045(1) 


0.4 


-0.00244(4) 


0.0030(1) 


0.5 


-0.00191(1) 


0.0022(6) 


0.6 


-0.001606(6) 


0.00176(2) 


0.7 


-0.001397(3) 


0.00145(1) 


0.8 


-0.001241(1) 


0.00124(1) 


0.9 


-0.001107(1) 


0.001051(9) 


1.0 


-0.000979(1) 


0.000872(3) 


1.1 


-0.000849(2) 


0.000710(8) 


1.2 


-0.000706(3) 


0.00052(1) 


1.3 


-0.000543(4) 


0.00033(3) 


1.4 


-0.000335(7) 


0.00007(1) 


1.5 


-0.00005(1) 


-0.0002(1) 


1.6 


0.00034(1) 


-0.0004(1) 


1.7 


0.00093(6) 


-0.0021(5) 


1.8 


0.0020(1) 


-0.02(3) 



TABLE II. Numerical results for the 2-loop coefficients c^ 1 ' ^ = V] c^' ^ an< ^ c ^'^ = 



APPENDIX A: CONSTRUCTION OF VERTICES 



Here we will explain the basic idea of expanding the expressions 1 / vXXX, appearing in 
the overlap-Dirac operator, in powers of g , using a procedure introduced by Y. Kikukawa 
and A. Yamada [35]. We also provide the expression for V^(k 1: k 2 ) further below. 

In an integral representation, the combination -j== can be written as 

r°° dt 1 



VYtX 



r°° dt 1 



In fact, Eq. (Al) is valid for any operator X provided that X^X has no vanishing eigenvalues. 
We begin the desired expansion of the overlap-Dirac operator in powers of g , by setting 

X^X = X%X +Z (A2) 
O(g° ) 

The first term of Eq. (A2) corresponds to the inverse fermionic propagator, while Z leads 
to the vertices; for our 2-loop calculation we need to write Z up to O(g^) 

z = (XqXi + xlx ) + (x t x 2 + x{Xi + x\x G ) + (x^x 3 + x\x 2 + x\Xi + xlx ) 

" v ' V v ' V v ' 

o(gh) o(sg) o( 9 l) 

+ (X t X 4 + X\X 3 + X\X 2 + Xlx x + X\X Q ) +O(g 5 ) (A3) 

Using the above equations, we write the denominator on the r.h.s. of Eq. (Al) as 

t 2 + xtx ~ T 2 + xjx ( 1 ~ z t* + xjx + Z t* + xlx Z t 2 + xjx + (A4) 

From this point forward it is easier to work in momentum space since, taking the Fourier 
transform, the denominator becomes diagonal 

F.T.I ^-r—] = - 1 - (A5) 

L t 2 + Xj*o t 2 + u 2 (p) 1 ; 

(co 2 (p) defined in Eq. (38)). Combining Eqs. (Al) and (A4) we derive the Taylor expansion 
of . in momentum space 

VxJx 

J x _ [°° dt 2ir5(p' -p) P dt 1 , 1 

V^tX^.T. ' 7T t 2 + UJ 2 {pi) J-oo 7T t 2 + UJ 2 {p') {P,P> t 2 + U 2 ( P ) 

r°° dt f°° dk 1 rr . . 1 r,,i s 1 /ax 

^-^/-ccC^Tp^^ (A6) 



The integral over t can now be performed by closing the contour around the upper 
complex t-plane. Considering, as an example the third term on the r.h.s. of Eq. (A6), the 
result is 



(, ° dt 



1 



1 



J-oo 7T t 2 +UJ 2 (p') t 2 +LU 2 (k) t 2 +UJ 2 (p) Uj(p')uj(k)uj(p) [u{p') +LU (Jfe)] [u(k)+u(p)] [(jj(p)+u(j/)] 

Similarly we integrate all terms of Eq. (A6) over t and this leads to Eqs. (38)- (42) and to 
the expression for V A (k\, k 2 ) which is presented below 



V A \k 1 ,k 2 ) 



d 4 k 3 d 4 k A d 4 k 5 1 



(2vr)4 (2vr)4 (2tt) 4 12 ^ u(k pi ) + co(k P2 ) 



n 



^u}(k pi )u(k P2 )u(k p3 )uj(k P4 )(u}(kp 1 )/6 + u(k PB )/30) x 



\peS5 



[(X^k^k^xlih, k^X.ih, h)xl(k 5 , k 2 ) XG (k 2 ) 

+x 1 (k 1 , h)xl(k 3 , k A )x x {k A , hhiih^iks, k 2 ) 

+X 1 (ki, h)x\(k 3 , k A )xo(k A )x{(k A , k b )X 1 {k 5 , k 2 ) 
+X l (k l M)xl{kz)X l (kMx\(k A ,k b )X l {k 5l k 2 ) 

+xo(h)xl(k 1 , h)Xi(h, k A )x\(k A , h)X 1 (h, k 2 ) 

~fi( ^C^M^Pa)^^) + u(k P3 fj J x 

+Xi(fc l5 Mx^M^ife fc 4 )x£(^)X 1 (A;4, MxJftO^te, k 2 ) 
+X 1 (k 1 , k 3 )xl(k 3 )X 1 (k 3 , k^xlihjX^h, k 5 )x\(k 5 , k 2 ) X0 (k 2 ) 
+X 1 (k 1 , k 3 )x\(k 3 , k A )xo(h)xl(k A , k 5 ) Xo (h)xl(k 5 , k 2 ) X o{k 2 ) 
+xo(h)xl(k 1 , k 3 )xo(ks)xl(k 3 , k A ) Xo (k A )xl(k A , k^X^h, k 2 ) 
+ X o(ki)x\(k u k 3 ) X o(k 3 )xl(k 3 , k 4 )X 1 (k i , k^xUk^X^h,^) 
+ X0 (ki)xl(k 1 , k 3 ) X0 (k 3 )xl(k 3 , k 4 )X 1 (k il k 5 )Xl(k 5 , k 2 ) X0 (k 2 ) 

+Xo(ki)xl(k 1 ,k 3 )X 1 (k 3 ,k A )xl(k A )X 1 (k A ,k 5 )xl(k 5 ,k 2 ) Xo (k 2 ) 
+ X o(ki)Xl(k 1 , fc 3 )Xi(fc 3 , k A )X\(k A , h) X0 (k 5 )Xl(k 5 , k 2 )xo(k 2 ) 
uj 2 (k Pl )uj{kp 2 )uj(k P3 ) 



+ E 



\peS 5 



uj(ki)u(k 2 )uj(k 3 )uj(k A )u;(k 5 ) 



(u(k P2 )[u(k p J/2 + u;(k P3 )/6}+u;(k p ^^ x 
Xo(*i)*l(fci, k 3 ) X0 (k 3 )xl(k 3 , k A ) X o(k A )x\(k A , k 5 )xo(h)xl(k 5 , k 2 ) Xo (fo 



(A7) 



(S5: permutation group of the numbers 1 through 5) 
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